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Abstract 

We propose a linear model for the dynamics of amplitudes associated with formation of phase 
slip centers in a one-dimensional superconductor. The model is derived taking into account the 
fact that, during the formation of phase slip centers the wave number of the superconducting 
phase remains practically constant. The model captures various forms of amplitude modulations 
associated with PSCs in closed analytic forms. 
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One dimensional (ID) superconductor has its width smaller than the Ginzburg-Landau 
(G-L) coherence length C,{T) (average separation between electrons in a cooper pair) where 
T is the temperature of the system. Near the critical temperature Tq of the superconducting 
(SC) to normal (N) phase transition, the coherence length ^{T) is quite large and a finite 
width channel or wire of superconducting material can be treated as a ID system. Such 
a system has the pure SC state (characterized by having no Ohmic voltage drop along the 
wire despite having current through it) globally stable below a critical current density jc 
(G-L critical current) at T < Tq- At a current density j > jc, being induced by an apphed 
voltage across the length of the wire, the SC state coexists with the N state up to an upper 
limiting current j2- Beyond j2 the N state is globally stable. Within the limits jc < j < j2, 
the SC state is metastable and gives up locally to the N state by occasional formation of 
phase slip centers (PSCs) . PSCs are the points in space where the amplitude (A) 

of the SC order parameter becomes zero. By the formation of PSCs, the complex order 
parameter of the SC state loses (or adds) a turn along the length of the wire and moves 
on to a state characterized by a smaller (or larger) wave number (q) to generally lower the 
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free energy (F) pj, 12|, y, Many localized forms of the SC amplitude and corresponding 
electrochemical potential /i (N state order parameter in the present context) have mostly 
been numerically observed in this regime 

mum 

under various circumstances. For 



ID superconductors, although it is an old classic problem, new exact results which were 
unknown previously are constantly being explored [l^ in view of the vast applicational 
importance of the subject. Despite so much work being done, till date, there is no single 
relationship that can at least partially capture the function space of the localized amplitude 
variations and corresponding chemical potential profiles in a general way for the resistive 
state. The present work is focused at achieving this goal. 

In this letter we propose a model to study the amphtude modulations associated with PSC 
formation. We extend the notion of PSC formation, as envisaged by Little [3], to hypothesize 
that, the wave number of the SC phase during a PSC formation remain constant till the 
amplitude vanishes locally. This is the main physical observation which would immediately 
imply, the corresponding phenomenological G-L free energy is effectively function of a single 
variable (amplitude of the SC phase). The simple model, we propose here, will be based 
on this new observation of effective uni-dimensional G-L free energy. We will solve this 
model to explicitly get the closed analytical form for dynamic and static PSCs and other 
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amplitude modulations. The above mentioned fundamental hypothesis of the constancy 
of the wave number during local variations of the amplitude allows us to derive a range 
of amplitude modulated SC and corresponding N phase profiles on the same footing of a 
central relationship that relates these quantities. Such a unified derivation of a range of exact 
analytic results are remarkable in view of the prevailing notion of the analytic intractability 
of the full nonlinear model (CGL) for such solutions. This point gets justified when we 
provide analytic expressions for numerically observed amplitude modulations by Kramer et 
al [sl in 1977 which, to our knowledge, has not analytically been accounted for till date. 

The dynamics of a ID superconductor in the resistive regime is given by CGL equation 
Q as 

ijjt + ifitp = + (1 - IV'HV' 

j = Im{tp* V ^) - f^x- (1) 

Here, we consider, the simplest possible form of the time dependent CGL to demonstrate 
our novel method of exploring the function space of localized amplitude modulations. In 
what follows, the treatment would be general, disregarding details of the dependencies on 
the coefficients and that is why we have set most of the coefficients to unity in Eq.l. In 
this model, ip is the superconducting order parameter which is complex valued. The system 
being ID, x is distance along the wire from some arbitrary origin. Let us consider, the length 
of the wire is L with cross sectional area a. The /i is electrochemical potential which can be 
considered as the order parameter of the N phase within the scope of G-L phenomenology. 
The j is current (density) through the superconducting wire and the suffix t and x of ifj 
indicate of the partial derivatives of the SC order parameter with respect to them. 

Eq.l has two stationary solutions. (1) ip = 0, fi = —xj (for constant j) which is the 
normal state and (2) = Ae*"^^, = 1 — A^, j = A^q, which is the superconducting state 
when fi = 5|. The SC order parameter can be visualized as a spiral wound around the x- 
axis (along the wire). If there are N turns along the length L, there is a total phase difference 
(j) = 27cN along L. Thus, the wave number q of the SC phase is a measure of number of turns 
the system has on a given length L, since, q = 2ttN/ L. The expression for the corresponding 
G-L free energy, when the steady state solution is purely superconducting, is given by 

F = La[{q^-l)A' + ^]. (2) 
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The form of F clearly indicates, the SC states with smaller q values are energetically favored. 
In other words, the spiraling SC order parameter would tend to lose its turns along the length 
of the sample to go to a lower free energy state. The turn loss happens when the amplitude 
of the SC phase locally vanishes and this local vanishing of the amplitude associated with a 
phase change of 27r is what known as PSC formation. The qualitative idea of PSC formation 
was originally given by Littel Q]. Little had pointed out that the boundary conditions make 
ip single valued and any internal (not at boundary) change of phase has to be achieved by 
a part of the spiral crossing the central axis that means a local vanishing of the amplitude. 

We argue that, so long as the amplitude does not exactly become zero at a point, while 
forming the PSC, the wave number of the SC phase remains necessarily unaltered because 
number of turns along the length of the sample remains unaltered. The width of a spiral 
can always shrink locally keeping turns unaltered. There is no apparent geometric reason 
as to why it should not happen. Moreover, the above mentioned amplitude wave number 
relationship is a characteristic feature of only the stationary SC state and it does not ensure 
a coupling between the amplitude and the wave number while a PSC forms. Our argument 
necessarily means a separation of the dynamics of the wave number and amplitude of the 
SC phase. So, we extend Little's idea to put forward the hypothesis that, the phase relax- 
ation only happens after the phase change is in place somewhere locally along the wire and 
amplitude and phase dynamics are independent of each other during such a process. 

At this stage, refer to Fig.l where a schematic diagram of the free energy F has been 
plotted against amplitude of the SC phase for the states with wave numbers go (upper 
branch) and — 2'k/\ (lower branch) where A is the wavelength of the SC phase with 1 
less turn than that having a wave number go- During this process of formation of PSC, the 
system will move uphill (induced by thermal fluctuations or whatever) from P to Q (along 
the upper branch in Fig.l) where the amplitude actually becomes zero. Now, being at Q the 
SC phase can only lose one turn and can eventually take the lower branch to come down to 
the minimum at R corresponding to the wave number go — 27r/A. Interesting to note that, 
the assumption of constancy of the wave number essentially makes the free energy defined 
in one dimensional space and thus, there are no other paths available for the system to move 
to the lower branch while forming PSCs. 

Therefore, through out this transition, the wave number remains fixed either at go or, 
on average to go — 27r/L (while relaxing) which makes the system remain confined to the 
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free energy range dominated by the quadratic terms of the Eq.2. At the time of relaxation, 
since it is a small change in the wave number that has to be accommodated, we consider it 
not to perturb the free energy profile in Fig.l appreciably. Note that, the position of the 
free energy trough moves towards smaller A values besides the depth getting reduced as q 
approaches unity (i.e. q becomes larger). Thus, an allowed variation in q could make the 
system move to the quartic term dominated part of the free energy while amplitude gets 
reduced. But, the effective constancy of q does not allow for such opportunity. So, any 
form of amplitude modulations consistent to this excursion be reasonably found out from 
the linear part of the dynamics (Eq.l). 

The linearized system that we would consider from now on, as our model, for looking at 
amplitude modulations is 

ipt + ifiip = ipxx + (1 - ^o)^', (3) 

which might be complemented by the same expression of current density as in Eq.l, but, we 
do not require any explicit mention of this other equation in what follows. Note that, one 
can not simply drop nonlinear term from Eq.l to get the linearized from of the dynamics. 
Instead, one has to replace the by the constant which is the amplitude of the pure 
SC phase at the free energy minimum. This form of the model preserves all the steady state 
properties of Eq.l. Recall the form of the pure SC solution of Eq.l where the super-current 
density (measurable quantity) can give the wave number of the SC state as j = q{l — q'^) and 
that also fixes the value of the constant Aq in Eq.3 as a function of the current density. To 
find out the amplitude modulations from the proposed model we would look at the dynamics 
of the amplitude at a larger scale. Practically, amplitude variations happen over a much 
larger scale (proportional to (^(T)) fs^) and the intended results would show up at these 
larger scales. 

Let us expand the order parameter ip a.s ip = A{X, t)[iIj^ + etp^] where the amplitude of 
the order parameter A{X, r) is a function of slow space and time scales X and r respectively. 
Take the N phase order parameter fi to be 0(e) and that is because we consider the associated 
fi profile to an amplitude modulation would be at the same scales. The small (< 1) parameter 
e is the ratio of the small to the large length scales i.e. the scales at which the phase of the 
stationary SC state and its amplitude vary respectively. Since, the dynamics of the system 
in this resistive state is slow, the only time scale under effective consideration here is the 
slow one t = er and we introduce multiple length scales as a; = a; + eX. The standard tools 
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of singular perturbations can now be brought to bear. Expanding up to 0(e) we get 



dr 



= A{X, t) 
+ e2A{X,T) 



dr 

dXdx 



+ e 



dA{X, r) 
dX 



dx 



dx 



+ {l-A^)A{X,T)[^'' + e^']. 



Thus, at 0(1), 



(4) 



(5) 



C + (1 - ^')^° = 0, 

is a zero eigenvalue equation and the solution is of the form ip^ = y4(X, r)e*'^*^"^''^^^. In the 
next order (0(e)), 

V^L + (1 - A')^' = i^^^ + r) 



,^ ,dA{X,T) . ,dq{X,T) 
- t2q{X,T) \' ' ^ixAiX.T) ' ' 



dX 



dr 



+ {2xq{X, r) - 2t)A{X, r) ^'^[^' )e''}ix,r). 

oX 



(6) 



The solvability, at this order, requires the coefficient of the secular term on the right hand 
side of the Eq.6 vanish. Considering the evolution of the wave vector to be independent of 
the amplitude dynamics we separate the equations for the amplitude and the wave vector 
as 



dA . , dA 

— + lefiA - i2q— = 

ot ox 

IX— + 2[xq - I)— = 0. 
ot Ox 



(7) 
(8) 



In Eq.7 and 8 the smaller scales have been restored which effectively introduces an e factor 
to the middle term of the Eq.7. 

The fixed points of the Eq.8 are q = qo (constant) and q = i/x. Considering q = qo(= 
1 — Aq), writing fi = —Jnx one can get a solution 15|] for the full time dependent form of 
the Eq.7 as 

A{x, t) = e"™^' [e^'^e^^* + e-^'^e"^^*] (9) 

The constants in the above dynamic solution m = ejN/Mo, ^ = ±2goA' and A' is a real 
constant. This solution explicitly represents a dynamic PSC of the system. The solution is 
real at the origin, where due to superposition, the oscillatory part becomes cos cut. Every- 
where else, for nonzero A', its a superposition of spirals (in time). A zero A' is not allowed 



for dynamic PSCs because uj is proportional to A' and no dynamic PSCs will form if A does 
not oscillate at the origin. Its important to note that, the oscillation frequency being pro- 
portional to wave number, the PSCs will be infrequent as the wave number approaches zero 
i.e. the system approaches the global minimum. This is not surprising, since, the height of 
the barrier it has to overcome is the same as the depth of the well d = ^(1 — ^q) and, given 
that, no saddle is available to bypass the barrier in ID. In a recent paper [l6|, dynamic PSC 
generation has been shown as a consequence of the PT (parity and time reversal) invariance 
of the CGL which is also evident from the form of the above mentioned PSC solution. The 
solution is invariant under the joint symmetry operations of parity and complex conjugation 
reflecting the necessity of having PT invariance in the CGL for dynamic PSCs. The temporal 
oscillation of the solution irrespective of explicit consideration of the thermal fluctuations is 
interesting to note. Such localized oscillations were noticed numerically by Kramer et al lal 
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and latter on (also numerically) by others (see [3|, |8|, |9| and references therein) within the 
resistive regime. These intrinsic dynamic PSCs are the dominant contributors to the phase 
loss at a region away from the Tc where large thermal excitations are less probable jsj. 

To look at the static exact solutions captured by our model, consider the time independent 
part of Eq.7. In a simplified form this static part of the Eq.7 



^JnA(.)^±-,(.), (10) 

relates the amplitude of the SC phase to the chemical potential and we argue that this is 
the central constraint that gives a selection for /i corresponding to a static modulation in A 
and vice versa. Following are the three solution pairs corresponding to some localized forms 
of the amplitude A 

1. A{x) = Aoe~"o^',/i = -^x 

2. A{x) = Aosechax,/! = — ^^tanhax 

3. A{x) = Aoe~"-°'^^ cosh ax, /i = ^ (— 2aoX + atanhax) 

In the first pair of conjugate forms of A and fi, a Gaussian localized bubble of the SC 
phase is present on a predominantly normal phase over the entire wire which carries a 
constant Ohmic current. These Gaussian bubbles of the SC phases are similar in form 
to .he _al„ found ou. enseal nuole, Wa«.ToMn a, fl. Suoh bubb,. of 
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SC phase are found to appear when one moves down the current from j2 towards jc and 
the pre-existing N phase loses stabihty to an SC phase. If the bubbles are bigger than a 
critical size (critical nucleus), they grow and push the N phase out towards the end of the 
wire, otherwise, these SC bubbles vanish and the N phase remains stable. In the second 
pair of solutions, the hyperbolic tangent profile of /i represents an intervening normal phase 
connecting two states characterized by different constant electrochemical potentials. Thus, 
a localized amplitude modulated SC phase is again shown to coexist with a normal phase 
that basically connects two reservoirs of charges at two different electrochemical potentials. 
In Fig. 2, we have plotted the third pair of the above solutions on the same scales as in the 
Fig.l in ref.jsl for a comparison. The figure in the ref. 5] has been obtained numerically 
by integrating a nonlinear model similar to Eq.l. However, to our knowledge, there is no 
analytical expression existing for such structures till date. The remarkable resemblance of 
the analytical plot to the numerically obtained structure supports our claim that Eq.lO is 
the constraint that the amplitude and the phase dynamics remain decoupled during such 
structure formations. The parameters that have been used here to plot the analytic solution 
are Aq = 0.6, 2go/e = 0.35, qq = 0.3 and a = 0.7. 

A static PSC can have the form of A(x) = Antanh Ixl, which corresponds to the form of 
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solution proposed by Langer and Ambegaokar for a small current density j [1|, y] . At the 
origin, where the PSC is sitting, corresponding fi{x) = ^[1/ tanhx — tanhx] will have a 
singularity. The analytic form we have derived here for the /i compares quite well with the 
numerically obtained one in ref. 31] (figure 18 in ref. 31]). Since, in the resistive regime, the SC 
and the N phase coexist, there could be some fronts present that separate these two phases 
jsl. For an example, consider an SC phase of the form A{x) = Aoe^/(l + e^). The Eq.3 
immediately produces corresponding solution for fi as fi{x) = ^(1 — 6^^/(1 + e"^')). Thus, an 
SC phase on the right hand side is separated by a corresponding N phase on the left and 
vice versa by these fronts. The central constraint (Eq.lO) correctly captures the numerically 
obtained structures in exact mathematical forms. 

To understand the implications of the existence of above mentioned amplitude modula- 
tions on the stability of the SC phase in the resistive regime, consider perturbing the fixed 
points of Eq.7 and 8, which are the above mentioned stationary local structures, with in- 
finitesimal perturbations 6 A and 6q. One can readily linearize Eq.7 and 8, in a general way. 
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for all of the above mentioned local structures to get 

dt \6q) \ ^Kl+^'-f)) \h 

where k is the characteristic wave number of the perturbations. The growth rates of the 
perturbations are Ai,2 = -(ze/io+^), («^-^)- Real part of the growth rate = is 
positive for opposite signs of k and go- What interesting about this linear stability analysis 
is that, it shows in the presence of these amplitude modulations the SC state is locally 
unstable to perturbations. The constant amplitude SC phase is known to be locally stable 
and globally unstable in the resistive regime. The allowed amplitude modulations facilitate 
transitions of the SC state to a lower free energy regime by making it locally unstable. 

To conclude, we have proposed a new model to find out exact analytic forms of the 
amplitude modulations associated with PSC formation and others. The model has been 
obtained on the basis of a hypothesis that the wave number of the SC phase remains constant 
during such structure formations. The solutions of the present model when compared with 
the already existing numerical results and experimental facts about localized modulations in 
the resistive phase vindicates our assumption of constancy of the wave number and resulting 
decoupling of the amplitude and phase dynamics. Linearity of the present model is an 
important aspect and its this linearity that allows us to get exact analytic solutions for a 
number of modulated states. We also have eventually derived a relation between the SC 
and N order parameters in the resistive regime which, we argue, is the basic constraint 
capable of producing all the observed stationary conjugate forms of the competing orders. 
An important consequence of the present analysis is that, it provides an exact measure of 
the free energy barrier in terms of the measurable super current density that the system is 
supposed to overcome while PSC forms. A revisit of the existing thermodynamics of such 
systems in the light of this new finding is worth doing in view of the large departures of the 
previously theoretically obtained numbers from the experimental ones. 
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A 

FIG. 1: Schematic diagram of the free energy F against amphtyde A. During the formation of a 
PSC when the system loses a turn to lower the free energy the path the system would follow will 
be from a to o along the upper branch and them come down to b from o along the lower branch of 
the diagram. 




11 



